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1. Introduction 

Spectroscopy with heavy quark is one of the major topics in high energy physics. Since the 
b quark is very heavy (~ 4.6 GeV), most theoretical approaches regard the b quark as a static, 
non-relativistic particle, or treat it with heavy quark effective theory or the relativistic heavy quark 
formalism. In other words, the heavy b quark is not treated as the Dirac fermion. This introduces 
complicated normalization procedures, in which the systematic errors are difficult to control. 

In QCD, all quarks are excitations of Dirac fermion fields, thus it is vital to preserve this 
important feature in any theoretical study. However, so far, it remains a challenge for a single 
lattice to accommodate such a wide range of quark masses from m u u quark to mt,. If we restrcit 
our study to the hadrons containing only s, c, and b valence quarks (ranging from 140 MeV to 4.6 
GeV) satisfying the constraints m q a < 1 and M^L > 4, then it is possible to accommodate these 
hadrons on a 32 3 x 60 lattice, with the inverse of lattice spacing a~ l ~ 7.68 GeV. 

In this proceeding, we report the first study of heavy meson spectra [|l|], by treating all heavy 
and light quarks as Dirac fermions, in lattice QCD with exact chiral symmetry. Namely, we restrict 
ourselves to mesons with quark contents bb, cb, sb, and cc. Our results of masses and decay 
constants of the pseudoscalar mesons B s and B c , and the masses of the vector mesons B* and B*, 
have been presented in [^]. 

We adopt the optimal domain-wall fermion (ODWF) proposed by Chiu [||, ^ as our fermion 
scheme, which possesses optimal chiral symmetry for any finite N s (the number of sites in the 5-th 
dimension), and also preserves the salient features of the Dirac fermion in the continuum. With 100 
gauge configurations generated with the single plaquette action at j3 = 7.2, we compute the point- 
to-point quark propagators for 33 quark masses in the range 0.01 < m q a < 0.85. We set the fifth 
dimension of ODWF, N s , to 128, and the stopping criteria of the conjugate gradient (CG) iterations 
to 10 , such that the chiral symmetry breaking a and the residual of the quark propagators e 
satisfy 

o = \Y t S 2 Y\/\Y^Y\ < KT 14 , £ = \\{D c + m q )Y - I|| <2x 10" 11 

where S is the optimal rational approximation of the sign function of 75D VV , D w is the Wilson Dirac 
operator plus a negative parameter —mo- Here (D c + m q )~ l is the valence quark propagator, and Y 
is the solution vector of solving the quark propagators via the CG algorithm. Note that our massive 
Dirac operator is exponentially local for all the quark masses. In Fig. [I], we plot the magnitude of the 
4-D effective Dirac operator of ODWF, D e ^(x;0) (i.e., the overlap-Dirac operator with Zolotarev 
approximation of the sign function) [||] along the t-axis, for one gauge configuration. Evidently, it 
is exponentially-local. Furthermore, we have checked that for any element of D ac ^^ (x; 0) (where 
a,b = 1,2,3 and a,j5 = 1,2,3,4), its magnitude is also exponentially-local along any direction, for 
all gauge configurations and all 33 quark masses. 
Next, we measure the time-correlation function 

C r (0 = (^Hr(D c + m Q )-lr(D c + m q )^}\ (1.1) 

for scalar (S), pseudoscalar (P ), vector (V), axial-vector (A), and tensor (T) mesons, with Dirac ma- 
trix r= {1, 75, Yi, 75 Yi, 75 74 Yi} respectively, for both symmetric (rag = m q ) and asymmetric (fixed 
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Figure 1: (a) The magnitude of the an element of D e ff(0,?;0,0) (with color index =2, and Dirac index = 1) 
from the origin to any point along f-axis, for one configuration, at m q a = m^a = 0.68. (b) Similar to (a), but 
averaging over all elements with different color and Dirac indices. 



niQ with various m q ) quark masses. Finally we extract their masses and pseudoscalar decay con- 
stants, compare our results with the experimental data, and also make some theoretical predictions. 

2. Determination of a , m c , m s , and nif, 

In Ref. [Q], we determine the inverse lattice spacing from the pion decay constant, with the 
experimental input f n = 131 MeV. However, in this work, our smallest quark mass turns out to 
be too heavy (about m s /2), thus the chiral extrapolation to m q ~ does not seem to be feasible. 
Instead, we use the mass and decay constant of the pseudoscalar meson T7 C (2980) to determine the 
bare mass of charm quark m c and cr 1 simultaneously, which can be seen in the following. 

We measure the time-corelation function Eq.(|l.l|) of the pseudoscalar meson Cp(t) (T = 75) 
for symmetric quark mass mq = m q , and fit it to the usual formula G(t), to extract the meson mass 
mpa and the decay constant fpa: 

G(t) = ^-[e-^ + e-" 1 "^}, f P a = 2m q a^ (2.1) 

Thus, we obtain the ratio of mp/fp for each m q a we have computed, which can be compared with 
the ratio of mass and decay constant of T] c . 

Although only the m^ c is measured in the high energy experiments, we can obtain the theo- 
retical value of niriJf^ = 6.8(2) from our previous study of pseudoscalar mesons on the 20 3 x 40 
lattice at j3 = 6.1 [Q]. We find that the closest value of mp/fp = 6.8(1) occurs at m q a = 0.16, 
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lattice 


MS 


MS 


PDG 


m q 


bare mass 


(at p = 2 GeV) 


(at ju = m q ) 


average 


m b 


5.22(4) GeV 


4.81(5) GeV 


4.65(5) GeV 


4.20(7) GeV 


m c 


1.23(4) GeV 


1.13(4) GeV 


1.16(4) GeV 


1.25(9) GeV 



Table 1: The bare mass (m q a ) of b and c quarks, and their masses in MS scheme at scales \i = 2 GeV 
and jj. — m q , respectively, in comparison with the PDG average in the last column. 



thus we fix m c a = 0.16. Further, using the experimental input m^ c = 2980 MeV, we determine the 
inverse of the lattice spacing a -1 = 7680(59) MeV. To check the goodness of our determination, 
we measure the time-correlation function of the vector meson c%c, and extract its mass equal to 
3091(11) MeV, which is in good agreement with 7/^(3097). 

Next, we determine the bare masses of strange and bottom quark, m s a and m^a, by extracting 
the mass of vector meson from the time-correlation function Cy(?) for various m q a we have com- 
puted. At m q a = 0.02, we obtain my = 1027(38) MeV, in good agreement with 0(1020); while at 
m q a = 0.68, we obtain my = 9453(3) MeV, in good agreement with T(9460). Thus we determine 
m s a = 0.02 and m^a = 0.68, respectively. 

Note that although the spatial size of our lattice (L ~ 0.8 fm) seems to be small, the finite size 
effects should be well under control since our pseudoscalar meson masses satisfies mpL > 4 even 
for our smallest bare quark mass m q a = 0.01. 

In order to compare our results of and m c with the high energy phenomenology, we have to 
obtain the lattice renormalization constants Z m = Zj 1 and transcribe them to the usual MS scheme, 
where Z s is the renormalization constant for yy. In general, the renormalization constants should 
be determined nonperturbatively. However, in this work, the lattice spacing is rather small (a ~ 
0.026 fm), we suspect that the one-loop perturbation formula [§] 

Z,(/x) = l + -f T [ln(flV) +0.17154], (mo = 1.30) (2.2) 

already provides a very good approximation for Z s . At j8 = 7.2, cT x = 7.680(59) GeV, our results 
for m^ s and m^ s at scales ji. = 2 GeV and their masses (/I = m q ) are listed in Table [l| Comparing 
with the PDG [Q] average, our result of m^ s (mb) turns out to be a little higher, but nJf s (m c ) is in 
good agreement with the PDG average. 



3. Charmonium cc and Bottomonium bb 



First of all, we check to what extent we can reproduce the charmonium mass spectra which 
have been measured precisely by high energy experiments. Our results of the mass spectra of the 
lowest-lying states of charmonium are summarized in Table ^[ The first column is the gamma 
matrix we use to compute the meson time-correlation function (Eq.(Ll). The second and the third 
columns are the J PC of the state and the conventional spectrascopic notation. The fourth column 
is the time range for fitting the time-correlation function to the formula in Eq.(2. 1). The # 2 /dof 
is listed in the fifth column, and the extracted meson mass is listed in the sixth column, where the 
first error is the statistical error, and the second error is the estimated systematic error based on all 
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r 


J PC 


n 2S+l Lj 


[tmimtmax] 


X 2 /dof 


Mass(MeV) 


PDG 


i 


0++ 


l 3 fl> 


[19,42] 


0.32 


3413(14)(9) 


X c o(3415) 


75 


o-+ 


1% 


[22,38] 


1.02 


2980(10)(12) 


t? c (2980) 


Yi 


1 


l 3 5i 


[19,38] 


1.18 


3091(11)(14) 


7/^(3097) 


75 Yi 


1++ 


l 3 A 


[18,42] 


0.51 


3516(13)(8) 


&i(3510) 


75747 


1+- 




[19,43] 


0.41 


3526(1 3X9) 


Ac(3524) 



Table 2: The mass spectra of the lowest-lying (n = 1) charmonium cTc states obtained in this work, in 
comparison with the PDG values in the last column. 



r 


J PC 


n 2S+l Lj 


[tmin j tmax] 


X 2 /dof 


Mass(MeV) 


PDG 


i 


0++ 


l 3 Po 


[21,39] 


0.21 


9863(15)(8) 


^o(9859) 


75 


o-+ 


1% 


[27,35] 


0.72 


9383(4)(2) 


T7fe(9300) ? 


Yi 


1 


1 3 S! 


[20,39] 


1.19 


9453(3)(2) 


T(9460) 


7s7 


1++ 


l 3 Pl 


[22,38] 


0.13 


9896(20)(8) 


Zm(9893) 


75747 


1+- 


l'A 


[22,38] 


0.10 


9916(22)(8) 





Table 3: The mass spectra of the lowest-lying (n = 1) bottomonium bTb states obtained in this work. The 
last column is the experimental state we have identified, and its PDG mass value. 

fittings satisfying # 2 /dof < 1.3 and |? max — t m i n \ > 6, with ? m j n > 10 and ? max < 50. The last column 
is the corresponding state in high energy experiments, with the PDG mass value. Evidently, our 
results are in good agreement with the PDG values. Note that our result of the hyperfine splitting 
(l 3 Si - 1%) is 111(14)(18) MeV, comparing with the PDG value 118 MeV. 

For the pseudoscalar T] c , we also obtain its decay constant f^ c together with its mass through 
Eq.([0]). Our result is fa = 438 ± 5 ± 6 MeV, where the first error is the statistical error, and the 
second error is the estimated systematic error, as described previously. 

Next, we turn to the bottomonium (bb) states. Our results of the mass spectra of its lowest- 
lying states are summarized in Table [| For the pseudoscalar rib, it was first reported by ALEPH 
Collaboration [gj], but it has not been confirmed by other HEP experimental groups. Our theoretical 
values turns to be a little deviate from their experimental result. In addition, we also determined its 
decay constant fa = 801 ±7 ± 5 MeV. 

Finally, we note that the tensor meson hb has not been observed in high energy experiments, 
thus our result of its mass serves as the first prediction from lattice QCD with exact chiral symmetry. 

4. Mesons with quark contents sb and cb 

The unitarity of CKM matrix is one of the crucial tests for the Standard Model. In order to 
extract the CKM matrix elements, the decay constants of heavy-light pseudoscalar mesons (e.g., 
/b, /b s , Id, and fu s ) have to be determined precisely. In Ref. [Q], we have presented our theoretical 
results of /o and fo s from lattice QCD with exact symmetry, which are in good agreement with the 
recent experimental results [Q] from CLEO Collaboration. 
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r 


J p 


n 2S+l Lj 


\tmin i tmax] 


Z 2 /dof 


Mass(MeV) 


PDG 


i 


0+ 


l 3 Po 


[20,40] 


0.38 


5852(15)(12) 


B; 7 (5850) 


75 


o- 


l l So 


[26,32] 


0.56 


5385(27)(17) 


£,(5368) 


Yi 


1- 




[25,33] 


0.50 


5424(28)(19) 


B*(5412) 


7s Yi 


1+ 


I 3 P { 


[18,42] 


0.68 


5884(16)(13) 




75 747 


1+ 


l l Pi 


[18,38] 


0.62 


5897(16)(12) 





Table 4: The mass spectra of the lowest-lying bTs meson states obtained in this work. The last column is 
the experimental state we have identified, and its PDG mass value. 
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J p 


n 2S+l Lj 


[tmin i tmax] 


X 2 /dof 


Mass(MeV) 


PDG 


i 


0+ 


l 3 Po 


[19,41] 


0.26 


6732(13X9) 




75 


o- 


I 1 So 


[19,38] 


1.14 


6278(6)(4) 


B c (6286) 


Yi 


1- 


1 3 S! 


[19,38] 


1.31 


6315(6)(5) 




YsYi 


1+ 


l 3 A 


[19,41] 


0.27 


6778(12)(7) 




YsYaY 


1+ 


l 1 A 


[18,42] 


0.45 


6796(10)(7) 





Table 5: The mass spectra of the lowest-lying bTc meson states obtained in this work. The last column is 
the experimental state we have identified, and its PDG mass value. 

In this work, we extract the masses and the decay constants of pseudoscalar mesons with sb 
and cb, via fitting their time-correlation functions Cp(t) to the formula 

F(t) = J—[e-^ + e- m ^ T -% fpa = {m +mQ ) a ^ (4 .1) 

Our results for B s are m Bs = 5385(27)(17) MeV, and f Bs = 253(8)(7) MeV. The former is in good 
agreement with the PDG value (5368 MeV). The later has not been measured in HEP experiments 
yet, thus our result serves as the first prediciton from lattice QCD with exact chiral symmetry. 

In Table |], we present our results of the mass spectra of the lowest-lying states of beauty 
mesons with quark content sb. Here we have identified the scalar bs meson with the state B* ; (5850) 
observed in HEP experiments, due to the proximity of their masses. Theoretically, this implies that 
6*^(5850) possesses J p = + , which can be verified by HEP experiments in the future. Further, 
our results of the masses of the axial-vector and tensor mesons, which have not been observed 
experimentally, serves a prediction from lattice QCD. 

For the pseudoscalar meson B c , we obtain m Bc = 6278 (6) (4) MeV, and f Bc = 489(4) (3) MeV. 
Our result of m Bc is in good agreement with the experimental value 6286(5) MeV measured by CDF 



Collaboration [jlOp, but f Bc has not been determined in HEP experiments. In principle, f Bc can be 
measured from the leptonic decay — > / + v/, since its decay width is proportional to f Bc \V c b\ 2 . 

In Table |5|, we summarize our results of the mass spectra of the lowest-lying states of mesons 
with beauty and charm. Except for the pseudoscalar B c , other states have not been observed in 
experiments. It will be interesting to see to what extent the experimental results would agree with 
our theoretical values. 
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QTq [t mi n,t max } X 2 /dof Mass(MeV) />(MeV) PPG 

by 5 b [27,35] 0.72 9383(4)(2) 801(7)(5) ^(9300) 

675c [19,38] 1.14 6278(6)(4) 489(4)(3) B c (6287) 

by 5 s [26,32] 0.56 5385(27)(17) 253(8)(7) 5,(5368) 

cy 5 c [22,38] 1.02 2980(10)(12) 438(5)(6) T} c (2980) 

Table 6: The decay constants of pseudoscalar mesons obtained in this work, together with their masses. 
They are identified with the corresponding PDG mesons listed in the last column. 

5. Concluding remarks 

We have performed the first study of treating b, c and s quarks as Dirac fermions in lattice QCD 
with exact chiral symmetry. The lowest-lying mass spectra of mesons with quark contents bb, cc, 
sb, and cb are determined, together with the pseudoscalar decay constants (which is summarized in 
Table ^). Furthermore, the b and c quark masses are determined (see Table [[]). Our results suggest 
that lattice QCD with exact chiral symmetry is a viable approach to study heavy quark physics 
from the first principles of QCD. 

For systems with u / d quarks, one may use several quark masses in the range m u j d <m q < m s 
to perform the chiral extrapolation. With a coarser and larger lattice, say, 42 x 64 at j3 = 7.0, it 
is possible to accomodate a wider range of quark masses m s /4 < m q < nib. This study is now in 
progress. 
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